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An exact non-perturbative calculation of the fourth-order 
anomalous correction to the scaling behavior of a random 
shell-model for passive scalars is presented. Importance of 
ultraviolet (UV) and infrared (IR) boundary conditions on 
the inertial scaling properties are determined. We find that 
anomalous behavior is given by the null-space of the inertial 
operator and we prove strong UV and IR independence of the 
anomalous exponent. A limiting case where diffusive behavior 
can influence inertial properties is also presented. 

Recently, there has been growing evidence, both nu- 
merically and experimentally, that fully developed tur- 
bulence is characterized by anomalous scaling of the ve- 
locity structure functions F p = (\v(x + r) — v(x)\ p ). In 
particular, it has been shown that, within the inertial 
range (i.e. 77 >C r <C L, 77 being the Kolmogorov dissi- 
pation scale and L the outer length of external forcing) 
F p (r) ~ r Zp , where z p is a non linear function of p which 



do not follow dimensional counting, i.e. z p ^ 



One 



OJ — _ F , q . 

of the most challenging scientific issue is to develop a 
theory which allows a systematic computation of z p by 
using the equation of motions. 

Recently [|l]-[| this issue has been addressed by study- 
ing a highly non trivial "toy model" introduced by 
Kraichnan, namely the advection of a passive scalar 
by a random, Gaussian velocity field, white-in-time 
and whose two-point velocity correlation function is 
given by:(v i (x,t)v j (x',t')) = 5(t - t')D v \\x - x'\), 
with Dij(x) — Dij(0) — Dij(x). Here, D is the d- 
dimensional velocity-field structure function : Dij (x) = 
Do\x\£ [(d — 1 + £)6ij — £xiXj\x\~ 2 ] , where the scaling 
exponent of the second order velocity structure function, 
£ (0 <C £ < 2), is a free parameter. Higher order velocity- 
field correlation functions are fixed by the Gaussian as- 
sumption. 

Although such a choice is far from being realistic, many 
interesting analytical and phenomenological results have 
been obtained for this toy-model. Due to the delta corre- 
lation in time, moment equations to all orders are closed. 
In (l]U, for the first time Kraichnan gave the closed ex- 
pression for S2 (r), where 



S p (r) = {\$(x)-0(x + r)\ p )~r c '. 



(1) 



6{x) being the scalar field transported by the turbulent 
velocity field. In the inertial range S p (r) ~ r^ p and the 
set of scaling exponents C p fully characterizes intermit- 
tency. In la theory for all structure functions is pro- 



posed and an explicit formula for £ p is derived. The main 
physical outcome is that all structure functions of order 
greater than two have intermittent corrections and that 
intermittency should be connected to some non trivial 
matching between advective and diffusive properties of 
the model. 

In y||] it has been shown that intermittency of the 
scalar structure functions is connected to the properties 
of the null-space of the linear operator appearing in the 
equation of multipoint passive scalar moments. More- 
over, in [|2|,|3| a perturbative expression of intermittency 
correction, as a function of the parameter £ and of the 
system dimensionality has been derived. Both in |jj and 
in [y,g| , matching conditions at infrared (IR) and ultravi- 
olet (UV) scales should be taken into account. The same 
problem for the case of passive vectors was addressed in 
. Finally, in all cases, universality in the scaling expo- 
nents is supposed to be preserved. 

In H two of us have worked out an even simpler toy- 
model which displays connections to the physics of a 
passive scalar advected by a random velocity field, be- 
ing at the same time more tractable both analytically 
and numerically. In [pi the intermittent properties of 
a shell-model for a passive scalar advected by a delta- 
correlated random velocity field have been investigated. 
The main results presented in [B were that: (i) the second 
order structure function has no anomalous scaling, i.e. 
C2 = 2 — £; (ii) all structure functions of order larger than 
two have anomalous corrections; (iii) anomalous behav- 
ior tends to vanish when approaching the laminar regime, 
£ = 2. For all the three points the model agrees with the 
Kraichnan model. 

The importance of UV and IR boundaries for the anoma- 
lous scaling was left unanswered in [g|. In particular, 
numerical simulations were unable to distinguish among 
contribution coming from the inertial null-space and pos- 
sible singular behavior introduced by the boundary con- 
ditions. 

In this letter we show how to compute exactly and non- 
perturbatively the inertial scaling behavior of the fourth- 
order structure functions. The main result is that the 
scaling properties are completely dominated by the null- 
space of the inertial linear operator and strongly uni- 
versal. The signature of UV and IR cut-offs is due to 
the presence of sub-dominant terms which weakly per- 
turb the pure-scaling behavior of the inertial operator. 



Anomalous scaling is calculated for any £ > 0. We also 
present some results which support the strong singular 
nature of the limit £ — * 0, In this limit, due to the non- 
local nature of interactions, it is not possible to neglect 
UV effects if the diffusive scale, kd, is taken fixed. 

We first introduce a simplified version of the shell 
model discussed in Ref. M . The model is defined in terms 
of a shell-discretization of the Fourier space in a set of 
wavenumber defined on a geometric progression k n = X n , 
with A > 1. Passive increments at scale r n = k~ x are de- 
scribed by a complex variable 9 n (t). The time evolution 
is obtained according to the following criteria: (i) the 
linear term is a purely diffusive term given by — Kk„6 n ; 
(ii) the advection term is a combination of the form 
k n 9 n 'U n "; (iii) interacting shells are restricted to nearest 
neighbors of n; (iv) in the absence of forcing and damp- 
ing the model preserves the volume in the phase-space 
and the passive-energy E = J2 n \9 n \ 2 . Properties (i), (ii) 
and and (iv) are valid also for the original equation of a 
passive scalar advected by a Navier-Stokes velocity field 
in the Fourier space, while property (iii) is an assumption 
of locality of interactions among modes. This assumption 
is rather well founded as long as <C £ < 2. The model 
is defined by the following equations (m = 1,2,...) 



<it 



K ^mWrn(t) 



i[c m 0* m+l (t)u* m {t) + b m e* m _ x {t)u* m _ x {t)) + S lm f(t) (2) 

where the star denotes complex conjugation and b m = 
—k m , c m = k m +i for imposing energy conservation in the 
zero diffusivity limit. Boundary conditions are defined as: 
Mo = #o = 0. The forcing term <5i m /(£)is gaussian and 
delta correlated: lf(t)f(t')) = FiS(t - t') acts only on 
the first shell. In numerical implementations, the model 
is truncated to a finite number of shells N with the ad- 
ditional boundary conditions 9n+i = 0. 
Following Kraichnan M we assumed that the velocity 
variables u m (t) and the forcing term /(£) are indepen- 
dent complex Gaussian and white-in-time, with scaling 
law: (um{t)u* n {t')) = S(t - t')8 nm d m , d m = k m ^ Due to 
the delta-correlation in time, we can close the equations 
of motion for all structure functions. Numerical simula- 
tions show that the model has the same qualitative in- 
termittency of the model studied in B . 

In this letter we concentrate on the non-perturbative 
analytic calculation of the fourth-order structure function 
Pmm = {{9m9* n ) 2 ) ^ ^m° 4 ( tne lowest order with non- 
trivial anomalous scaling). The closed equation satisfied 

by P mg = {{o m o; n ){e q e* q )), is, 

Pmq = {5l,m Em +$l,qE q )Fi - n(k m + k 2 q )Pmq 
+ [-PmqC 2 m d m ((1 +S q . m+1 ) + X^ 2 (l + Sg.rn-i)) 

\*m+l,qCrn™''>T*>\'®q-,'>'n) * *~Tn— l,q^rn ^"m-~ 1 \^ * Oq,m) 

+(q «-► m)] (3) 



where E n 
(ft as: 



j#*). We can symbolically represent eq. 



P — T 

1 mq — -^"1 



mq,lp*lp 



kV 



mq,lp*lp 



Pin + T, 



mq 



(4) 



where X and T> are the inertial and the diffusive 4-order 
tensor and J- is the forcing term. 

Our main result is derived by using the following ansatz: 
the symmetric matrix P mq , which fully determines the 
scaling properties for any fourth-order quantity in the 
model, can be described as: 

P n ,„+i = CiP n m {I > 0), P n ,„_i = DlP n , n (I > 0) (5) 

The independency of C; and Di from n is equivalent to 
demand absence of strong boundary effects, i.e. the ma- 
trix is formally infinite-dimensional. Clearly this must 
be verified a posteriori showing that the solution we are 
going to present is UV and IR stable. 
Using (0) we obtain: 



C, 



i+i 



D 



i+i 



Qd D l D 1 

which is equivalent to write: P n +i,n+i = k l P n ,n where 
Ci/Di — fc; and £4 = 2(2 — £) — p±. As usual we indicate 
by p4 the anomalous correction to the scaling exponent. 

Let us notice that (||) does not force the solution to 
have global scaling invariance: only the diagonal part is 
requested to have pure scaling. 

Let us proceed by analyzing (g) restricted to the in- 
ertial operator and for the diagonal (m = q) and sub- 
diagonal terms (q = m — 1): 



Pn 



IP 



,c- m d m {-l-x + 2{C 1 + D 1 x)) 



*m,m—l — ^m 



X 



m,Tn—l^rn' m '\ 

c 2 



x z Ci 
R 



c\ 



(6) 



(7) 



where we have posed x — A^~ 2 and R = C\/D\. By 
plugging the scaling (||) in (0) one obtains two equation 
in three unknowns which can be taken to be C\ and the 
ratios CijC\ and R. Numerical investigation suggest the 
following " scaling ansatz" : 



P n ,n+l = CiP n m, with C; = Cikf_? 



and 



Pn,n-l = DiP n , n , With Di = Dxk^l 2) Pi 



(8) 



(9) 



where anomalous correction is felt only in the IR part 
of the matrix. By plugging this scaling in (M) we end 
up with two eqs. in two unknowns and we can calcu- 
late pi. Let us anticipate that this (wrong) assumption 
gives results in very good agreement with the numeri- 
cal simulations, indicating that the true solution is not 



very far from having pure scaling behavior. In order to 
solve the full problem, without imposing any "pure scal- 
ing" behavior, we analyze the other entry of the ma- 
trix P nq with q ^ n and q =/= n — 1. Let us put 
ll =Di +1 /D h 5 l =C l+1 /Ci. 

It is then possible to show that for I > 1 by plugging 
the scaling (H) in the inertial part of (0) and studying the 
equation for P„ ]Ti ±i we obtain two recursion equations: 

- (1 + a;)(l + x l ) + 7/ (i? + x l+1 ) + — (x l + ~) = (10) 

7/-i ti 

(1 + x){\ + x l ) - \^6i(R + x l+1 ) - ^(x l + ~) = (11) 

These two relations can indeed be seen as two maps con- 
necting successive values of 7; and Si respectively. By it- 
erating forward (backward) the map (not) we move from 
the diagonal (IR boundary) to the IR boundary (diago- 
nal) along a row of the matrix Pi >n . By iterating forward 
(backward) the map (O) we move from the diagonal (UV 
boundary) to the UV boundary (diagonal) along a row 
of the matrix P/. n . 

Let us first note that the two maps are not indepen- 
dent, i.e. they satisfy our scaling ansatz Si = R^i and 
therefore we are going to consider only one of the two 
in what follows. In order to test stability under weak 
perturbation of boundary conditions in the map ( |l0| ) we 
are interested in the behavior by backward iterations, i.e. 
iterating from I = oo to I = 0. In the limit / — > oo and 
£ ^ 2 the map (|10| ) has only two fixed point correspond- 
ing to 7^ = x/R and 7! = 1/R. It turns out that 7^ 
is stable for back iterations, i.e. iterating from the IR 
boundary (I 3> 1) to the diagonal (1 = 0). The global so- 
lution can now be obtained by a self-consistent method. 
First, let us take as initial value for R the value that 
one would have guessed from imposing " pure scaling" as 
discussed previously, then we can iterate ( |l0|) from the 
boundary toward the diagonal and finding the value for 
C2/C1. This value can be used to close (0) exactly. Next, 
with the improved value for R, one can restart the full 
procedure getting a new improved value of R and so on 
up to the moment when the new value of R reach it fixed 
point. 

In figure 1 we show the computation of p± obtained by nu- 
merical integration of equation (0) as a function of £. In 
the same figure we plot the P4 as a function of £ obtained 
by the analytical solution previously discussed. As one 
can see, the agreement is perfect. Let us notice that it 
is impossible to go by numerical simulations to values of 
£ very near zero because of strong diffusive effects which 
completely destroy scaling behavior. 
Let us notice that it is the strong stability under UV and 
IR perturbations that allows us to iterate consistently 
the procedure. We have therefore proved that anomalous 
scaling comes only from the inertial operator and that it 
shows a very strong degree of universality as a function 



of the forcing and dissipative mechanisms, at least as 
far as the situation with £ = const. > and molecular 
diffusivity n — > is considered. 

Some new and interesting phenomenum happens when 
we are in the other possible asymptotic limit (£ — > 0) 
and fixed molecular diffusivity. A simple dimensional 
argument tell us that the following relation holds: 



const. k. 



(12) 



where with kd we mean a preassigned diffusive scale such 
that all the inertial dynamics is at k <C kd- One can show 
that in such a situation the diffusive operator V in (Q) 
gives a contribution of the form 



*n,n — const. (fc n / kd) k n ± n ,n 



(13) 



Therefore, the diffusive perturbation is absolutely neg- 
ligible in the case when £ =/= is fixed and kd — * 00 
while it becomes a singular perturbation when kd is kept 
fixed and £ — * 0. A detailed analysis of relations (Q) 
shows that in the singular limit £ — > there appear an 
infinitesimal interval of values of £ ~ where p± —* 0, i.e. 
the anomalous correction tends to vanish with a partic- 
ular shape which depends on the constants appearing in 
relation (|12|). 

This concludes the presentation of our results. 
Let us remark that the perfect agreement of our inertial 
null-space solution with the numerical simulation per- 
formed with finite diffusivity and in presence of forcing 
is the clear demonstration that the scaling behavior is 
completely dominated by the inertial operator. For any 
finite system, IR and UV effect weakly perturb the pure 
scaling solution. Our result shows that the inertial opera- 
tor is perfectly suitable for picking all anomalous aspects 
but in the case where strong non local interactions (dy- 
namically produced) completely destroy inertial proper- 
ties introducing diffusive effects at all scales (£ — > and 
kd fixed). 

Discussions with B. Dubrulle, G. Falkovich, U. Frisch, 
I. Procaccia and M. Vergassola are kindly acknowledged. 

Figure caption 

Figl: analytical ansatz (continous line) and numerical 
results (squares) for ^4 are plotted for various values of 
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